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A method and application results are described for using mathematical models in statistical analysis and con¬ 
trol of the annealing regime of a glass band transported inside a tunnel furnace. The results of a computational 
experiment estimating the efficiency of the algorithms of the annealing regime control are discussed. The pos¬ 
sibility for reducing waste and improving the quality of polished glass annealing is demonstrated. 


The introduction of the ISO-9002 quality control system 
in the production of sheet glass implies the use of statistical 
methods in controlling the process of polished glass produc¬ 
tion. We will consider the methods and results of the applica¬ 
tion of mathematical models for analysis and control of the 
annealing conditions of a moving glass band in a tunnel 
furnace. 

The method consists of the following stages: 

- analysis of the annealing process with the aim of 
choosing an impulse for evaluation of the quality of anneal¬ 
ing and refining the determining factors; 

- development of “regime - quality” mathematical mo¬ 
dels describing the dependence of the glass quality parame¬ 
ters on the annealing temperature regime and the monitored 
variable input parameters; 

- development of an algorithm for adjusting the tempera¬ 
ture regime depending on the geometrical dimensions of the 
band and the line capacity, taking into account the controlled 
disturbances. 

The performance of technological line 1 of an annealing 
furnace produced by the KNUD company was analyzed. The 
line capacity is 420tons/day, the glass band width up to 
3300 mm, and the band thickness up to 6 mm. For one month 
the annealing regime data were registered every 4 h and the 
glass quality parameters were determined. 

The quality of annealing was estimated based on the re¬ 
sidual stress value, the flat band curvature, and the percent of 
glass waste after annealing. 

The residual stresses in glass were monitored at five dif¬ 
ferent points across the band width. The measurement data 
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are correlated, and the correlation coefficient is 0.83 - 0.95. 
This made it possible to estimate the residual stress based on 
the measurement of stresses in the middle of the band. To im¬ 
prove the accuracy, the residual stress was calculated as the 
arithmetic mean of five measurement results. 

The curvature of the manufactured glass was monitored 
in samples taken from the left and from the right sides of the 
glass band. The lateral and longitudinal curvatures of the 
lower and the upper band surfaces were measured. The mea¬ 
surement results correlate, and the correlation coefficient 
ranges from 0.60 to 0.89. The generalized parameter of the 
band curvature after annealing was calculated as the arithme¬ 
tic mean of 8 sample measurements. 

The daily glass waste after annealing was calculated 
based on the operating time of the glass-crushing machine 
considering its efficiency. 

To refine the list of developed models, the closeness of 
the correlation between the indicators of the produced glass 
quality was determined. The matrix of pairwise correlation 
coefficients is given in Table 1. 

The data in Table 1 point to a significant correlation 
existing between the parameters of glass annealing quality. 
The negative correlation between the glass waste and the re¬ 
sidual stress is evidence of a substantial amount of waste 


TABLE 1 


Glass annealing 
parameters 

Residual 

stress 

Curvature 

Glass waste 

Residual stress 

i 

-0.48 

-0.32 

Curvature 

-0.48 

1 

0.25 

Glass waste 

-0.32 

0.25 

1 
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Furnace length, m 

Fig. 1. Averaged curves of annealing temperature variations at dif¬ 
ferent sections across the glass band: 1 ) left edge; 2 ) center; 3 ) right 
edge. 

generated in producing thin glass, and the positive correla¬ 
tion between the waste and the curvature indicates an in¬ 
crease in the glass waste resulting from high-curvature 
glasses. The negative correlation coefficient between the cur¬ 
vature and the residual stress reflects the greater curvature in 
thin glass, which have lower residual stresses. The signs of 
the correlation coefficients correctly characterize the rela¬ 
tionships between the individual parameters of glass anneal¬ 
ing quality. As the coefficients of correlation between indi¬ 
vidual quality parameters are low in modulus, they served as 
the basis for developing all models describing the quality of 
annealing. 

Depending on the geometrical dimensions and the che¬ 
mical compositions of glass, a specific optimum annealing 
regime is selected for each type of glass. The present study 
considers the case where the chemical composition of glass 
is relatively stable, and its variations are neglected. 

The annealing regime depends on the temperature varia¬ 
tions in the tunnel furnace and the velocity of the band mo¬ 
tion. The temperature inside the annealing furnace is moni¬ 
tored in zones A, B, C, and D along the tunnel. The tempera¬ 
ture variations in each zone are registered at five points 
across the band width. To choose an impulse characterizing 
the annealing temperature conditions, the temperature varia¬ 
tions along the furnace axis measured at various points 
across the band were subjected to statistical analysis (Fig. 1). 

A confidence interval for the mean temperature in dif¬ 
ferent zones was constructed for each curve. The averaged 
curves of temperature variations along the annealing furnace 
at different cross-sections across the glass band fall within 
the overlapping confidence intervals. This makes it possible 
to describe the variations in the mean temperature by a linear 
dependence of the type 

Q = a 0 -a l x, (1) 

where a 0 is the free coefficient of the equation; a x is the coef¬ 
ficient of the straight line slope; x is the coordinate of the po- 


TABLE 2 


Parameter 

Residual stress 

Curvature 

Glass waste 

a o 

0.03 

0.13 

0.09 

a l 

-0.27 

0.07 

-0.02 


sition of the thermocouple in different furnace zones along 
the longitudinal axis. 

Formally the free coefficient corresponds to the glass 
band temperature at the entrance to the annealing furnace 
(x = 0), which should depend on the temperature of the free 
end of the glass band at the exit from the float tank. The 
closeness of the correlation between a 0 and the temperature 
of the exit end of the band is estimated by a correlation coef¬ 
ficient equal to 0.07, which makes it possible to neglect this 
dependence. The coefficient of the straight line slope charac¬ 
terizes the intensity of temperature variation along the longi¬ 
tudinal axis of the annealing furnace, which determines the 
process of annealing of the moving glass band. The effect of 
a 0 and a 1 on the quality of annealing is characterized by the 
correlation coefficients shown in Table 2. 

As can be seen, a significant correlation exists only be¬ 
tween a 0 and the glass curvature, as well as between a 1 and 
residual stress. The effect of the coefficients on other para¬ 
meters can be neglected due to its low significance. 

To select the factors affecting glass annealing, the fol¬ 
lowing monitored variables that correlate with the glass qua¬ 
lity parameters were considered: the temperature of tin in the 
float tank at span 15, the temperature of the glass band end at 
the exit from the float tank, the glass density, the production 
line efficiency, and the thickness and output rate of the pro¬ 
duced glass band. Due to the analytic dependence of the effi¬ 
ciency on the band parameters and output rate, the efficiency 
was excluded from the number of analyzed factors. The glass 
thickness and the output rate were left for analysis. The tem¬ 
perature of the exit end of the band significantly depends on 
(correlates with) the temperature of the tin at span 15. The 
pair correlation coefficient for these temperatures is 0.72, 
which makes it possible to consider only one temperature, 
i.e., the temperature of the end of the band at the exit from 
the float tank. 

Thus, five factors were selected for the construction of 
regression models: the temperature of the glass band end at 
the exit from the float tank, the glass density, the band thick¬ 
ness, the output rate, and the parameters characterizing 
the temperature distribution along the annealing furnace 
( a 0 and a x ). 

In developing regression models of the annealing pro¬ 
cess, the conditions of glass melting and the glass band for¬ 
mation over the tin melt, as well as the chemical composition 
of the glass, were assumed to be constant. These assumptions 
simplify the construction of models but affect their precision. 
The justifiability of these assumptions was substantiated by 
subsequent studies. 

The regression models were constructed from a sampling 
of 180 experiments reflecting the annealing furnace perfor¬ 
mance during one month. The structure of the models was re¬ 
fined using the procedure of consecutive regression analysis. 
At first, all the above listed variables were included in the 
original structure, and then the insignificant variables were 
consecutively eliminated using the Student /-criterion. The 
elimination continued, while the multiple correlation coeffi- 
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cient kept growing and the residual dispersion of the model 
kept decreasing. As a result, the following regression models 
with estimates were obtained: 

- residual stresses in glass (pm/mm): 

a = b 0 + 5 + b 2 v - b 5 a l , R 2 = 0.86; 

- the glass band curvature (mm/mm): 

C=-b 1 5 + b 3 Q exit -b 4 p, R 2 = 0.87; 
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Fig. 2. Conditions of glass annealing on LPS-1 line: 1 ) actual va¬ 
lues; 2 ) model values. 


- glass waste (%): 

W = - b x 8 + b 3 9 exit - b 4 p + b 5 aj, R 2 = 0.74, 

where b 0 - b 5 are the regression coefficients; 5 is the band 
thickness, mm; v is the band velocity, m/min; R 2 is the mul¬ 
tiple coefficient of linear correlation; a l is the parameter of 
model (1) of annealing temperature, K/m; 9 exit is the tempe¬ 
rature of the glass band end at the exit from the float tank, 
°C; p is the glass density, g/cm 3 . 

Let us consider the efficiency of using the developed 
models for statistical analysis and control of the glass band 
annealing process. The matrix of connectivity of the exit pa¬ 
rameters of annealing with the annealing temperature condi¬ 
tions, the produced band parameters, and the glass density is 
shown in Table 3. 

Among the specified factors, the intensity of temperature 
variation along the axis of the annealing furnace was chosen 
as the controlling factor. As the chilling intensity grows, the 
waste at the exit from the furnace decreases. At the same 
time, the residual stress in the glass increases. The glass band 
curvature is virtually not affected by variations in the chilling 
intensity. The other factors represent monitored variables and 
cannot be controlled in the course of annealing. Thus, the 
glass band thickness, the output rate, and the band tempera¬ 
ture at the exit from the float tank are determined by the 
glass-molding technology and depend on the geometrical di¬ 
mensions of the band and the efficiency of the production 
line. The glass density depends on the batch composition and 
the glass melting process. Thus, one can only control the an¬ 
nealing process through varying the intensity of chilling the 
moving glass band. 

Considering the above, the annealing process control cri¬ 
terion can be written as a penalty function: 

F2 = X t | min (0.25 - Amal, 0) | + X 2 I min (a adm - c, 0) | + 
k 3 | min (0.025 - C, 0) | + k 4 W, (2) 

Oadm = 5+2.5(S-2), 

where Lj - A, 4 are the penalty coefficients; Amal is the daily 
variation in the intensity of glass band chilling in the anneal¬ 
ing furnace; cr adm and a are, respectively, the admissible and 
the actual residual stress in the glass band. 


TABLE 3 


Glass 

annealing 

parameters 


Determining factors 


band 

thickness 

band 

velocity 

chilling 

intensity 

temperature 
of the band 
exit end 

glass 

density 

Glass 

waste 

Inverse 

depen¬ 

dence 


Linear dependence 

Inverse 

depen¬ 

dence 

Residual 

stress 

Linear dependence 

Inverse 

depen¬ 

dence 



Curvature 

Inverse 

depen¬ 

dence 



Linear 

depen¬ 

dence 

Inverse 

depen¬ 

dence 


The first summand in the penalty function (2) limits the 
daily variation in the chilling intensity of the glass band in¬ 
side the annealing furnace, and the second one limits the 
value of the admissible residual stress in the glass, which de¬ 
pends on the produced glass thickness. The third summand is 
related to the glass band curvature, which should not exceed 
0.025. The final summand is proportional to the glass waste 
value. 

The penalty coefficients were selected experimentally: 
kj = 100, L, = 3, k 3 = 1, k 4 = 5. Simulating modeling of the 
algorithm of the glass annealing control was carried out us¬ 
ing the selected penalty coefficients on the LPS-1 production 
line at the Borskii Glass Factory between January 1 and De¬ 
cember 31, 2000. 

The annealing temperature regime can be controlled us¬ 
ing the residual stress value. To accomplish this, the penalty 
coefficients A, 3 and X 4 in expression (2) are taken equal to 
zero. In adjusting the annealing regime, the temperature of 
the exit end of the glass band did not vary and was equal to 
the actual temperature. The intensity of the temperature vari¬ 
ation along the annealing furnace axis is shown in Fig. 2. It 
can be seen that in statistical controlling, the temperature in¬ 
side the annealing furnace in producing 3-5 mm glass 
should be maintained virtually invariable. In the case of ma¬ 
nual control, a lower intensity of glass chilling intensity was 
maintained. 

The simulation modeling of the algorithm of statistical 
control of residual stresses in glass revealed the possibility of 
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Fig. 3. Temperature in the annealing furnace: 1 and 2 ) actual and 
model temperature in zone A; 3 and 4 ) actual and model temperature 
in zone B; 5 and 6) actual and model temperature in zone C; 7 and 
8 ) actual and model temperature in zone D. 

decreasing the average annual value of residual stresses in 
glass from 8.98 to 8.38 pm/mm, which is equal to 6.7%. In 
this case the mean quadratic deviation of residual stresses in 
produced glass remains virtually constant: 2.31 and 
2.06 pm/mm. 

The problem of controlling the annealing furnace opera¬ 
tion can also be stated as the selection of a regime ensuring 
the smallest amount of glass waste while satisfying the 
boundaries imposed on the admissible residual stress value 
and the glass band curvature. In this case, penalty function (2) 
is used. 

The variations in the annealing temperature are shown in 
Fig. 3. The model temperature in zone A virtually coincides 
with the temperature registered in the manual control of the 
annealing process, whereas the temperature in zone B is 
slightly lower. The temperatures in zones C and D differ in¬ 
significantly from the temperatures registered in the manual 
handling of the process. Only in the second half of the year 
were these temperatures lower. 

The chilling intensity in the model annealing regime dif¬ 
fers to some extent (by 9.3%) from the manual control re¬ 
gime. A faster chilling was observed in the second half of the 
year. The chilling regime in modeling is more stable, and the 


TABLE 4 



Actual residual stresses 

Model residual stresses 

Date 

mean 

monthly 

value 

mean 

quadrati 

deviation 

mean 

monthly 

value 

mean 

quadratic 

deviation 

23.01.00 

10.66 

1.68 

10.09 

1.60 

23.02.00 

7.99 

1.91 

7.14 

2.00 

24.03.00 

9.22 

1.75 

8.77 

1.34 

23.04.00 

9.47 

2.19 

7.79 

1.60 

23.05.00 

9.57 

2.69 

8.89 

2.44 

22.06.00 

7.76 

1.74 

7.34 

1.72 

21.07.00 

8.89 

1.88 

6.67 

1.13 

21.08.00 

8.65 

2.29 

8.16 

1.57 

20.09.00 

8.63 

2.33 

8.81 

1.70 

20.10.00 

8.08 

2.26 

8.3 

1.95 

19.11.00 

8.76 

1.6 

8.88 

1.51 

19.12.00 

8.87 

2.63 

9.14 

2.13 

Mean annual 
value 

8.88 

2.08 

8.33 

1.73 


mean quadratic deviation of chilling intensity is by 9% 
lower. 

The comparative results of this control algorithm model¬ 
ing and the manual control are shown in Table 4. 

The use of the annealing control algorithm makes it pos¬ 
sible to diminish residual stresses in glass after annealing by 
6%, whereas the spread in residual stress values decreases by 
16.8%. The curvature of the produced glass band virtually 
does not change in transition from the manual annealing con¬ 
trol to annealing using the control algorithm. The mean cur¬ 
vature is equal to 0.0150 and 0.0156 mm/mm, respectively. 

The average daily duration of the operation of the 
glass-crushing machine in the manually controlled annealing 
process was 3.45% of the day duration with the mean qua¬ 
dratic deviation equal to 2.7%. In glass annealing using the 
control algorithm, the average daily operation of the glass¬ 
crushing machine was shortened to 2.69%. This yielded an 
increase by 0.78% in the glass output on the LPS-1 line, 
which is equivalent to an additional 1169.9 tons of glass 
(three days of LPS-1 line operation). 

The performed studies confirmed the possibility of re¬ 
ducing the glass waste generated on the LPS-1 line and im¬ 
proving the glass quality by using the developed mathemati¬ 
cal models and control algorithms. 














